Cooper problem in the vicinity of Anderson transition 
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We study numerically the ground state properties of the Cooper problem in the three-dimensional 
Anderson model. It is shown that attractive interaction creates localized pairs in the metallic 
nonintcracting phase. This localization is destroyed at sufficiently weak disorder. The phase diagram 
for the derealization transition in the presence of disorder and interaction is determined. 
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The pioneering experimental results for normal-state 
resistivity of high temperature superconductors demon- 
strated a striking correlation between the optimal doping 
with maximal T c in the superconducting phase and the 
Anderson metal-insulator transition (MIT) in the nor- 
mal phase obtained in a strong pulsed magnetic field 
0. More recent experiments on the superconductor- 
insulator transition (SIT) in three dimensions (3D) 
which were done in various materials at different dopings 
and magnetic fields, also reveal close correlation between 
these transitions even if it is possible that the normal 
state remains metallic in some materials These ex- 

perimental results put forward the important theoretical 
problem of interaction effects in the vicinity of Ander- 
son transition in 3D. However, the full understanding of 
this problem is very difficult since even the origin of the 
high-T c phase is not yet established completely. Due to 
that it would be interesting to understand the effects of 
interaction and disorder in a more simple model of gener- 
alized Cooper problem || of two quasiparticles above the 
frozen Fermi sea which interact via the attractive Hub- 
bard interaction in the presence of disorder. In spite of 
apparent simplicity of this problem it is rather nontrivial. 
Indeed, even if the great progress has been reached re- 
cently in the investigation of localized one-particle eigen- 
state properties , the analytical expressions for the in- 
teraction induced matrix elements in the localized phase 
and in the MIT vicinity are still absent. Furthermore, the 
recent results for the problem of two interacting particles 
(TIP) in the localized phase demonstrated that the inter- 
action effects for excited states can qualitatively change 
the eigenstate structure leading to the appearance of de- 
localization ||7|— |To|l . Due to that the investigation of the 
ground state properties of the above model in the vicinity 
of the Anderson transition in 3D represents an interest- 
ing unsolved problem which can shed light on the origin 
of SIT in the presence of disorder. 

To investigate the above problem we study numerically 
the ground state properties of two particles with Hubbard 
on site attraction (U < 0) in 3D Anderson model at 
half filling. In this case the one particle eigenstates are 
determined by the Schrodinger equation 




E n 1p n + V(lp a -l + 1pn+l) = Eip r 



(1) 



FIG. 1. Probability distributions for TIP in the ground 
state projected on (x, j/)-plane: one particle probability f p 
for interaction U = — 4V (left column), interparticle distance 
probability f p d for U = —4V (centrum column), f p for [7 = 
(right column); the disorder strength is W/W c = 1.1 (upper 
line), W/Wc = 0.5 (middle line), W/W c = 0.3 (bottom line). 
All data are given for the same realisation of disorder for 
the system size L = 16 (see text for details). Upper line 
corresponds to the insulating noninteracting phase while two 
others are in the metallic one at U = 0. 

where n is the site index on the 3D lattice with periodic 
boundary conditions applied, V is the nearest neighbour 
hopping and the random on-site energies Ei are homo- 
geneously distributed in the interval [—W/2,W/2]. It 
is well known that at half filling (the band center with 
E = 0) the MIT takes place at W c /V » 16.5 with the 
insulating and metallic phases at W > W c and W < W c 
respectively (see e.g. jll],[l2]] ) . To study this problem with 
interaction it is convenient to write its Hamiltonian in the 
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basis of noninteracting eigenstates of the Anderson model 
that gives 



m, ,m 



^ mi , 7?i2, m 1 ,m 2 / V-m 1 ,m 2 



(2) 



Here Xmi ,m 2 are eigenfunctions of the TIP problem writ- 
ten in one-particle eigenbasis </> m with eigenenergies E m . 
The transition matrix elements Q m2 m > m > are ob- 
tained by rewriting the Hubbard interaction in the nonin- 
teracting eigenbasis of model ([!]). The Fermi sea is intro- 
duced by restricting the sum in (0) to m 1 2 > with un- 
perturbed energies E > > Ep. The value of the Fermi 

m l,2 

energy Ep ~ is determined by the filling factor fi which 
is fixed at fj, = 1/2. To have more close similarity with 
the Cooper problem we also introduce the high energy 
cut-off defined by the condition 1 < m 1 +m 2 < M. Such 
a rule gives an effective phonon frequency wd k M/L 3 
where L is the linear lattice size. Since the frequency u>d 
should be independent of L we keep the ratio a = L 3 /M 
constant when varying L. The majority of data are ob- 
tained for a w 30 [|l3) but we checked that its variation 
by few times did not affect the results. Due to on-site 
nature of the Hubbard interaction only symmetric con- 
figurations are considered. 

In fact the first studies of the model (Q) with the frozen 
Fermi sea had been done by Imry || with the aim to an- 
alyze the derealization effect of TIP in the proximity of 
Fermi level at finite particle density. This model was also 
studied numerically in where it was shown that near 
Ep the interaction becomes effectively stronger compar- 
ing to the ergodic estimate used in However the 
above studies ||Jl4|] were concentrated on the properties 
of excited states in the repulsive case U > 0. On the 
contrary here we analyze the ground state properties for 
the attractive case. Since U < 0, then even in the limit of 
large system size L the particles are always close to each 
other in the ground state that is qualitatively different 
from the case U > 0. In this way the model (Q) repre- 
sents the generalized Cooper problem in the presence of 
disorder. 

To study the characteristics of the ground state \m\ ,m 2 
we diagonalize numerically the Hamiltonian (|2J) and 
rewrite the eigenfunction in the original lattice basis |n) 
with the help of relation between lattice basis and one 
particle eigenstates |n) = ^ m Rn,m4>m- As the result 
of this procedure we determine the two particle proba- 
bility distribution F (111,112) in the ground state (here 
111,2 mark the positions of the two particles), from which 
the one particle probability /(ni) = X)n 2 ^( n i> n 2) 
and the probability of interparticle distance /d(r) = 
F(r + n 2 , n 2 ) with r = ni — n 2 are extracted. For 
graphical presentation these probabilities are projected 
on (x,y)-plane that gives f p (n x ,n y ) = £\ f{n Xl n yi n z ) 



and f P d(x,y) respectively. The typical examples of pro- 
jected probability distributions f p and f p d for different 
values of disorder W are shown in Fig. 1 . They clearly 
show that in the presence of interaction the ground state 
remains localized not only in the noninteracting local- 
ized phase (W > W c ) but also in the phase delocalized 
at U = (W < W c ). However the localized interacting 
phase abruptly disappears if disorder W becomes smaller 
than some critical value W S (U) < W c . For W < W s the 
ground state becomes delocalized over the whole lattice. 
At the same time the peaked structure of the interparti- 
cle distance distribution f p d clearly shows that the par- 
ticle dynamics remains correlated. In this sense we can 
say that the pairs exist for any strength of disorder but 
for W > W s they are localized while for W < W s they 
become delocalized. We assume that such a transition 
should correspond to the transition from insulating to 
superconducting phase in the many-body problem. 
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FIG. 2. Dependence of IPR £ on the rescaled disorder 
strength W/W c for U = (o), U = -2V (□), U = -W (o), 
U = -61/ (A) (empty/full symbols are for L = 10/i = 12). 
Insert shows the data for U — —41/ at L — 8 (x), 10 (o), 
12 (full diamond), 14(*); £ obtained in the Cooper approx- 
imation (see text) for U = —41/, L = 14 is shown by (+). 
Statistical error-bars are smaller than symbol size. Lines are 
drawn to adapt an eye. 

To analyse this transition in a more quantitative way 
we determine the inverse participating ratio (IPR) £ for 
one particle probability : l/£ = (X) n / 2 ( n ))' where 
brackets mark the averaging over 100 disorder realisa- 
tions. Physically, £ counts the number of sites occupied 
by one particle in the ground state. Its variation with 
system size L is shown in Fig. 2 for different strength of 
interaction and disorder. This figure shows that in the 
localized interacting phase W > W S (U) the £ value re- 



2 



mains finite and independent on size L while in the delo- 
calized phase it grows proportionally to the total number 
of sites L 3 . To find the critical disorder strength W s we 
compare the relative change of £ with L (8 < L < 14) 
with its relative change for the noninteracting case at 
the critical point W = W c . Then W S (U) is defined as 
such a disorder at which the relative variation of £ at 
\U\ > becomes larger than in the case U = 0. We 
note that near the transition the change of £ with W is 
so sharp that the derealization border is not really sen- 
sitive to the choice of definition. We also checked that 
the change of uj£> does not affect significantly the border 
W S (U) fig] . The phase diagram for SIT defined in the 
way described above is presented in Fig. 3. It shows that 
the interaction makes localization stronger so that the 
localized phase penetrates in the noninteracting metallic 
phase. However for sufficiently weak disorder derealiza- 
tion takes over. Qualitatively we can say that the attrac- 
tion creates a pair with a total mass (m p ) twice larger 
than the one particle mass and due to that the critical 
disorder strength becomes twice smaller (W s /W c ~ 0.5) 
since the effective hopping V e ff oc \/m p . Of course this 
argument is not sufficient to explain the exact border 
W S (U) obtained numerically but it gives a reasonable es- 
timate in the case of strong interaction. Further studies 
are required to explain the form of the border. 
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FIG. 3. Phase diagram for transition from localized (in- 
sulating I) to delocalized (superconducting S) phase in the 
ground state of the generalized Cooper problem (^|). Verti- 
cal dashed line shows the Anderson transition in absence of 
interaction. 

Another interesting physical characteristic is the cou- 
pling energy A of two particles in the presence of interac- 
tion. Its value is equal to A = 2Ep — E g where E g is the 



ground state energy in the presence of interaction and 
2Ep is equal to E g at U — 0. In the standard Cooper 
problem A > is related to the BCS gap and determines 
the correlation length of the pair. It is interesting to un- 
derstand how A varies with the disorder strength W at 
fixed interaction U. This dependence is presented in Fig. 
4. It clearly shows that A grows significantly with the 
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FIG. 4. Variation of the ground state coupling energy A 
in the model ([]) with the rescaled disorder strength W/W c 
for L — 12 and different interaction U = — 2V (full box), 
U = -W (full diamond), U = -6V (full triangle); the open 
symbols are for the same U values but in the Cooper ap- 
proximation (see text). The vertical dashed line W/W c = 1 
marks the MIT at U — 0; the other vertical dashed line 
W/Wc ~ 0.35 marks approximately the SIT line from Fig. 
3 for \U\ > 2. Insert shows the case U = -4V for L = 10 
(o), L = 12 (full diamond), L — 14 (*); dotted line shows the 
Cooper approximation case with U = —4V and L = 12 from 
the main figure; logarithm is decimal. 

increase of W at constant interaction U. We attribute 
the physical origin of this growth to the fact that at 
stronger disorder the rate of separation between particles 
becomes smaller that enhances enormously the interac- 
tion between them, hence A, as it was discussed in [ pT| . 
The dependence of A on W is changed drastically near 
W s that is related to the derealization transition. 

On the same figure we compare the exact value of 
A, found numerically in the model (^|), with its value 
Ac obtained by the Cooper approximation (mean field 
value). In this approximation only the matrix elements 

Qmi,m 2 ,m[,m' 2 with TOl = m 2 and TO 'l = m 2 are ke Pt in @ 
that corresponds to the original Cooper ansatz ||. The 
comparison shows that at weak disorder Ac is very close 
to exact A (see insert where dotted line coincides with full 
diamonds for W/W c < 0.35) while when approaching the 
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Anderson transition and beyond it (W/W c > 0.35) Ac 
becomes much smaller than A. This leads to the conclu- 
sion that the nondiagonal matrix elements (mi ^ mi and 
m 1 =/= rn 2 ), neglected in the Cooper approximation, play 
an important role near MIT. This is also clear from Fig.l 
according to which the localized states exist in the non- 
interacting metallic phase while according to the Cooper 
approximation pairs should be delocalized for W < W c . 
Indeed, for example the graphical image as in Fig. 1 
shows that for U = -4V, W = 0.5W C > W s the probabil- 
ity f p obtained in the Cooper approximation from (Q) is 
completely delocalized contrary to the real case in which 
the ground state is localized (Fig. 1). In addition to 
this case the average IPR within the Cooper approxima- 
tion is much larger than its real value obtained without 
approximation (see insert in Fig. 2). 

The data presented in the inserts of Fig. 2 and 3 clearly 
show that in the localized interacting phase W > W s 
the lattice size is sufficiently larger than the localiza- 
tion length and the values of £ and A correspond to the 
limit L — > oo |l6}| . At the same time for weak disorder 
W/W c < 0.2 the asymptotic value of A is very small 
and very large values of L are required to reach it. Such 
large L are also desirable to see better the propagation of 
pairs with large size. The main increase of A takes place 
in the metallic noninteracting phase at W s < W < W c . 
However in this region the TIP pair remains localized 
due to interaction that does not allow to obtain a gain 
in the value of A ~ T c . It would be interesting to find 
some possibility to dclocalize the pair in this region and 
to keep large A at the same time. 

In conclusion, our numerical studies of the generalized 
Cooper problem in the presence of disorder show that in 
the ground state the attractive interaction leads to lo- 
calization of pairs inside noninteracting metallic phase, 
contrary to the Cooper ansatz. This localization how- 
ever disappears at sufficiently weak disorder. The phase 
diagram for the transition to delocalized states is deter- 
mined as a function of disorder and interaction. 

We thank V.V.Flambaum, K.Frahm and O.P.Sushkov 
for stimulating discussions, and the IDRIS in Orsay and 
the CICT in Toulouse for access to their supercomputers. 



363 (1998) ( |cond-mat/9807332|). 

K. Segawa and Y. Ando, |cond-mat/9908124 



[4] 

[5] L.N. Coope r, Phys. Rev. 104 1189 (1956). 
[6] A. Mirlin, |;ond-mat/9907126| (to appear in Phys. Re- 
ports) 

[7] D.L. Shepelyansky, Phys. Rev. Lett. 73, 2607 (1994); 
p.201 in 0. 

[8] Y. Imry, Europhys. Lett. 30, 405 (1995); p.211 in 0. 
[9] D. Weinmann, A. Miiller-Groeling, J.-L. Pichard, and 

K. Frahm, Phys. Rev. Lett. 75, 1598 (1995); p. 221 in 

@- 

[10] Correlated fermions and transport in mesoscopic systems, 
Eds. T. Martin, G.Montambaux and J.Tran Thanh Van, 
Editions Frontieres, Gif-sur-Yvette, (1996). 

[11] B. I. Shklovskii, B. Shapiro, B. R. Sears, P. Lambrianides 
and H. B. Shore, Phys. Rev. B 47, 11487 (1993). 

[12] T. Ohtsuki, K. Slevin, T. Kawarabayashi, Ann. Physik 
(Leipzig) 8, 655 (1999). 

[13] In this case we found that lu d /V » 0.5W/W c + 0.1. 

[14] Ph. Jacquod and D.L. Shepelyansky, Phys. Rev. Lett. 
78, 4986 (1997). 

[15] For example W 3 /W c ~ 0.35/0.36 for a = 30/15 at U = 
-41/. 

[16] We checked that for W > W s , e.g. W = 12V, U = -4V, 
the excitated states with E — 2Ef < also remain local- 
ized, J. Lages and D.L. Shepelyansky (in preparation). 



http:/ /w3-physthco.ups-tlsc.fr/^dima| 



[1] G.S. Boebinger, Y. Ando, A. Passner, T. Kimura, M. 
Okuya, J. Shimoyama, K. Kishio, K. Tamasaku, N. 
Ichikawa, and S. Uchida, Phys. Rev. Lett. 77, 5417 
(1996). 

[2] V.F. Gantmakher, L.P. Kozeeva, A.N. Lavrov, D.A. 
Pushin, D.V. Shovku n, and G.E. Tsydy nzhapov JETP 
Lett. 65, 870 (1997) ( |cond-mat/9705197| ). 

[3] V.F. Gantmakher, M.V. Golubkov, V.T. Dolgopolov, 
G.E. Tsydynzhapov, and A. A. Shashkin, JETP Lett. 68, 



4 



